Ab Initio Investigation of Collective Charge Excitations in MgB2 
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A sharp collective charge excitation is predicted in MgB2 at ~ 2.5eV for q perpendicular to 
the boron layers, based on an all-electron analysis of the dynamical density response within time- 
dependent density functional theory. This novel excitation, consisting of coherent charge fluctuation 
between Mg and B sheets, induces an abrupt plasma edge in the experimentally-observable reflec- 
tivity. The existence of this mode reflects the unique electronic structure of MgB2, that is also 
responsible for strong electron-phonon coupling. By contrast, the acoustic plasmon, recently sug- 
gested to explain the high Tc, is not realized when realistic transition strengths are incorporated. 



PACS numbers: 74. 70. Ad, 78.20.Bh, 78.20.-e 



The high critical temperature, Tc ^ -iOK, of the newly 
discovered ^ conventional superconductor MgB2 has at- 
tracted much attention ||^-^ . A consistent picture within 
BCS theory pO|-p^ seems to be given by the qualitative 
analysis based on the large hole density of states at 
the Fermi surface, provided by the boron cr-bands. On 
the other hand, other possible contributions to pairing 
[f7|,p[p^ have not been totally ruled out, since a quanti- 
tative analysis of Tc, via realistic treatment of screened 
electron-electron and electron-phonon interactions [Q, 
has yet to be performed. The insight into the behavior of 
the cuprate superconductors resulting from consideration 
of the spin susceptibility suggests that it might be espe- 
cially interesting to examine the possibility of unusual 
collective charge excitations in MgB2 with such realistic 
treatments. In particular, a calculation of the dynamical 
density response would give detailed information on elec- 
tronic excitations, which control most of the macroscopic 
properties, including the dielectric function and thereby 
superconductivity p2| . 

Even within phonon- mediated theories, there currently 
exists large speculation of the value of the Coulomb pseu- 
dopotential n*, from 0.02 to 1 ||,^. An intriguing pro- 
posal Q has recently suggested that an anomalously 
small fj,* may originate from the existence of an acous- 
tic plasmon (AP) |l4|jr^ ], which develops from two char- 
acteristic charge carriers in MgB2: "heavy" holes of B 
CT-bands and "light" electrons and holes of B 7r-bands. 
This is an important issue, as the existence of an AP 
can drastically alter the dynamical electronic screening 
and may even provide an additional pairing channel . 
Thus, it is vital to re-examine the existing tight-binding 
result via a more realistic investigation that takes effects 
of the crystal potential fully into account. 

In this Letter we present results of an ab initio, all- 
electron, calculation of the dynamical density-response 
function of MgB2 within the framework of time- 
dependent density functional theory (TDDFT) Our 
key result is the prediction of a sharp collective mode 



whose origin is the strong coherent charge fluctuation be- 
tween parallel sheets of B and Mg. This novel charge ex- 
citation, found to reside within the [0002] zone-boundary 
gap, embodies a remarkable signature of the electronic 
structure of MgB2, and results in a dramatic change in 
the reflectivity near 2.5eV for q perpendicular to the 
boron layers. It thus explains naturally the strong in- 
crease in the Raman efficiency reported very recently . 
The impact of this mode on the Coulomb pseudopoten- 
tial will be briefly addressed. In addition, with the matrix 
elements that control the transition probability fully in- 
cluded, the recently proposed acoustic plasmon is not re- 
alized in our results, as a consequence of weak intraband 
transition of the "heavy carriers" . Finally, comments on 
the proposed two-band-type superconducting instabili- 
ties will be given to further illustrate the importance of 
including realistic transition strengths. 

Within TDDFT, the dynamical density response func- 
tion x{xjt; x' ,t') — 5p{x,t) / 5vcxt{x' ,t') with density p 
and external potential Vexti can be obtained through the 
following integral equation p6[ : 



(1) 



where x^^i^^^T^'^^') is the response function for 
"non-interacting" Kohn-Sham (KS) electrons, v{x — 
x') is the Coulomb interaction, and fxc{x,t; x' ,t') = 
5vxc{x,t) / 5p{x' ,t') accounts for dynamical exchange- 
correlation effects (where Vxc is the time-dependent 
exchange-correlation potential in TDDFT.) Working in 
Fourier space, can be calculated with the KS 

eigenenergies, ^ and the corresponding eigenstates, 

\k,n), of the ground state ||l7[|: 
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X {k, n\e-''^^+^>'^\k + q, n'){k + q, n'\e'^'i+'^'>^\k, n), (2) 
where G is a vector of the reciprocal lattice, n is a band 
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index, the wave vectors k and q are in the first Brillouin 
zone, and V is the normahzation volume. Equation (1) 
is then numerically solved as a matrix equation |Q . 

We stress that this formalism is rigorous. Even though 
a formal connection between the KS eigenenergies and 
the quasi-particle energies has not been established, the 
physical meaning of the KS band structure (including 
the empty states) is to be realized through its contribu- 
tion to the density fluctuation in ■ The only physical 
approximation introduced in this work is the functional 
form of Vxcixjt), for which the adiabatic local density 
approximation (ALDA) is chosen. Since only the 
long wavelength limit is of interest in the present work, 
we further drop fxc from Eq. (1), as its strength (within 
ALDA) cannot compete with the singular v without siz- 
able crystal local-field effects |Q. In this limit, the di- 
electric function e = {l + vx}^^ can be visualized directly 
from the single-particle transitions of the KS particles, 
with the relation e = 1 — vx^^ ■ 

Numerically, in order to resolve the features in the nar- 
row energy range where the AP may occur, the broad- 
ening parameter, commonly employed [T9| ] in place of 0^ 
in Eq. (2), needs to be very small(~ lOmeV). This in 
turn requires a dense fc-mesh, to carefully sample the 
momentum-energy phase space near the Fermi surface. 
The results presented in this work correspond to a uni- 
form mesh of 20 x 20 x 60. In addition, summing the nar- 
row (5-functions is avoided by first evaluating Eq. (2) on 
the imaginary frequency axis, followed by analytic con- 
tinuation via Pade approximants pO] . 
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FIG. 1. Calculated loss function (upper panel), dielectric 
function (middle panel), and reflectivity (lower panel), as well 
as the jellium counterparts (dashed lines) for q = 0.12A~^, 
along the c-axis. 

The calculated dielectric function and corresponding 
loss function (— Ime^^ — —vlmx) along the c-axis, 
which should be directly probed by inelastic scattering 
experiments of electrons or x-rays pi[| , are shown over a 
broad energy range in Fig. 1. The wide structure Q at 
uj ^ 20eV in the loss function is the conventional plas- 
mon corresponding to a density of ~ 8 electrons per unit 
cell {vs ~ 1.8ao). It suffices to note the correspondence 



between the near vanishing of the dielectric function and 
the peak line shape of the loss function, illustrating the 
collective nature of this over-damped mode. 

The most striking feature of the loss function shown 
in Fig. 1 is the sharp peak that occurs at ~ 2AbeV . 
The origin of this peak is evidently the prominent ab- 
sorption at beV in Ime, whose physics will be addressed 
shortly. Indeed, the strength of this charge fluctuation 
(notice the vertical scale of e) is large enough to gener- 
ate, through screening, an additional zero in Re e thereby 
inducing a new collective mode. The extremely long life- 
time (Aw < IQmeV) of this mode is a consequence of 
both the strong dynamical screening in the vicinity of its 
energy, as reflected in the large slope of Re e [0|3| , and 
the lack of particle-hole decay channels, courtesy of the 
large energy gap to be discussed below. 

This additional collective excitation has a dramatic 
physical impact on the reflectivity, as shown in the bot- 
tom panel of Fig. 1. When possessing enough energy to 
excite this mode, the incident light suddenly penetrates 
the surface of the material with no reflection, much more 
abruptly than the case of normal metal with a plasmon. 

The reduction of reflectivity also provides a natural 
qulitative explanation for the large enhancement of the 
Raman effeciency ||2^ in a recent 2AleV Raman scat- 
tering study of the E2g phonon, as more photons are 
allowed to enter the surface to interact with the system. 
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FIG. 2. Calculated band structure and density. The ar- 
rows in the band structure (left panel) indicate the large SeV 
gap between paraUel bands. The isosurfaces of density for 
states marked by the upper/lower arrows are shown in the 
upper/lower right panel, which contains 2x2x1 unit ceUs 
with Mg in the corners. 

Examination of the band structure, shown in the left 
panel of Fig. 2, reveals that the beV structure in Ime 
originates mainly from the interband transitions between 
parallel bands. These transitions are facilitated by the 
strong [0002] crystal potential, which reflects the layered 
structure of this system and gives a large beV gap (in- 
dicated by the arrows) at the F point. Intriguingly, as 
shown in the right panel of Fig. 2, where the real-space 
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charge distribution of states marked by the down and up 
arrows is illustrated, the states involved in the transitions 
reside in different planes and correspond to a B 7r-state 
and a layered conducting state of Mg s-symmetry. Evi- 
dently, the above-mentioned additional mode consists of 
coherent charge fluctuations between B and Mg sheets, 
under the strong dynamical screening of low-energy in- 
traband transitions. 

Note that this mode is of very different nature from 
the d-plasmon observed in a few post-transition metals. 
The d-plasmon consists of direction-insensitive collec- 
tive dipole-\ike oscillation of the localized electrons, orig- 
inating from transitions of flat d-bands to a wide range 
of final states. Thus, the corresponding absorption in 
Inie is usually 5-10 times smaller than the current case. 

The wave vector dependence of the new collective 
mode, controlled by the competition between dynamical 
intraband screening and interband charge fluctuation, is 
shown in Fig. 3, in which the surface illustrates underly- 
ing single-particle transitions given by Im e{q/ /[OOQl],uj). 
The first branch at ^ 0.5eV of the low-energy ridge of 
the surface corresponds to the intraband charge fluctua- 
tion of the heavy carriers (j)x,y cr-bands) to be discussed 
below, and the second branch to the light carrier (pz 
TT-bands), while the high-energy ridge demonstrates the 
beV interband transition. In the valley between these 
two main structures, where the phase space for single- 
particle excitation is closed due to the [0002] crystal po- 
tential, lies the weakly dispersed, induced collective mode 
(marked by dots), which eventually decays into electron- 
hole pairs of the light carriers at g O.GA"-*^. 




FIG. 3. Dispersion of collective modes and the sin- 
gle-particle decay channel. The dispersion of the in- 
duced/nominal collective mode is denoted by dots/circles un- 
til it becomes difficult to identify the peak energy of the plas- 
mon. The surface illustrates lme{q,u>) for (f//[0001]. 

It is natural to wonder about the impact of this addi- 
tional mode on the superconductivity, through reduction 
of the Coulomb pseudopotential, n*. Within the static 
approximation adopted by previous authors |l^,|l^, this 
new mode is found to have almost no effect on the static 



screened interaction, due to its very sharp line shape and 
small spectral weight compared to the conventional plas- 
mon. In fact, the static limit of our resulting screened 
interaction is almost identical to that obtained from the 
homogeneous electron gas, despite its highly anisotropic 
energy dependence. One thus would expect a fi* of "nar- 
mal" magnitude, if it is evaluated from the screened in- 
teraction However, due to the strong fc-dependence 
of the Fermi surfaces, even the very idea of treating the 
screening through a single Coulomb pseudopotential may 
not be adequate (c.f. Eq. (2.26) of Ref. 

Now we address the proposed ||] AP with calculated 
dielectric function for small momentum/energy, shown in 
Fig. 4. The two main peak structures in Ime are related 
to the unusual characteristic of the Fermi surface j|] . The 
peak of lower energy corresponds to the intraband charge 
fluctuation of the heavy carriers, while the main peak re- 
sults from the light carriers. More specifically, for a small 
momentum transfer q along the c-axis, the cylindrical 
Px,y tr-hole Fermi surfaces contribute to smaller en- 
ergy (heavier mass), compared to the three-dimensional 
Pz TT-hole and 7r-electron Fermi surfaces. This double- 
peak structure is precisely what led to the suggestion 
of AP 1^, as the plasmon energy of the heavy carriers 
might be greatly screened by the light carriers and re- 
duced to the upper edge of low-energy peak, whose en- 
ergy increases linearly with small q. 
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FIG. 4. Calculated dielectric function of MgB2, for small 
q's along the c-axis. The low/high-energy peaks of Im e corre- 
spond to charge fluctuation of the heavy/light carriers. Note 
the large difference in the spectral weight and that Ree 2> 1 
between these two main structures. 

However, this simple picture neglects electron-hole 
decay channels and the relative strength (or spectral 
weight) of these two carriers. As shown in Figure 4, 
for small q there is a large supply of charge fluctua- 
tions (Ime ^1) involving the light carriers at the upper 
edge of the low-energy peak. As a result, any collective 
mode of this energy would rapidly decay into electron- 
hole pairs (Landau damping) and have a negligible life- 
time. In fact, our analysis suggests that a gap in the 
light-carrier bands at the Fermi energy is needed to allow 
an AP to exist. Furthermore, in this energy range, Ree 
is considerably larger than unity, reflecting the fact that 
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the spectral weight from the heavy carriers is insignifi- 
cant compared to the strength of the screening from the 
light carriers. In short, the potential acoustic plasmon 
is totally screened out and should not contribute to /i*. 
(The low energy intra-cr-transitions may still lead to non- 
adiabatic corrections to the high energy, small q phonons, 
in particular the E2g modes, which are most important 
for superconductivity.) 

This kind of interplay between these two structures 
in e cannot be faithfully estimated by just counting the 
phase space, as performed in Ref. Q with parameter- 
ized tight -binding bands. Instead, a careful evaluation 
of the corresponding transition probability, controlled by 
the matrix elements in Eq. (2), is necessary. The com- 
monly made approximation of unit probability risks pro- 
moting unphysical charge fluctuations in the calculation, 
especially in the case of multi-band models . 

For a similar reason, the recently proposed two-band- 
type superconducting instability in MgB2 is also not 
supported by our results. For q ~ (0,0,7r/c), where the 
nesting between B 7r-electron and 7r-hole Fermi surfaces 
becomes "perfect" , there is no apparent enhancement at 
zero energy in our dielectric functions {q ~ 0.9A^^ in 
Fig. 3), contrary to the large interband polarizability ob- 
tained in Ref. 0. Even though the formalism adapted 
in this work is different from that in Ref. J? , based on 
the experience gathered in the past years 2^ , we are 
confident that our results should not suffer from qual- 
itative error for a simple sp system like MgB2. Thus, 
the enhancement of interband pair scattering in Ref. |^] 
is an artifact of omitting the transition probability (c.f. 
Eq. (9) of Ref. 0) as discussed above. 

In conclusion, a novel sharp collective mode, involving 
coherent charge fluctuation between B and Mg sheets, is 
predicted. This additional collective mode has a great 
impact on the optical properties of the material - sud- 
denly switching off the reflectivity in the corresponding 
energy. It thus naturally explains the strong increase 
in the reported Raman efficiency at 2AleV, by increas- 
ing chance of interaction. By contrast, the recently sug- 
gested AP is not supported by careful evaluation of the 
dynamical density response functions, with the full in- 
formation of all-electron KS band structure. The weak 
charge fluctuation from heavy carriers is totally screened 
out by the light carriers. Thus, AP is effectively ruled 
out as a possible mechanism that brings an anomalously 
small Coulomb pseudopotential, /i*. Our results further 
demonstrate the importance of inclusion of realistic crys- 
tal potential effects, through the band dispersions (gaps) 
and matrix elements, in the calculation of charge fluctu- 
ation for systems with non-spherical Fermi surfaces. 
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